defaultEncoding <- "UTF8"
library(dplyr)
library(sf)
library(terra)
library(raster)
library(DT)
library(ggplot2)
library(plotly)
library(leaflet)
casos <-
st_read(
"/vsicurl/https://marcelocubero.github.io/capas_proyecto/casos.geojson",
quiet = TRUE
)
casos <-
casos %>%
st_transform(4326)
casos_ <-
read.csv( file= "https://marcelocubero.github.io/capas_proyecto/casos_.csv",
sep = ",", encoding = "UTF-8")
fallecidos <- rast("/vsicurl/https://marcelocubero.github.io/capas_proyecto/fallecidos.tif")
fallecidos_rl <- raster::raster(fallecidos)
bins <- c(10, 100, 500, 1000, 4000, 7600)
pal <- colorBin("Blues", domain = casos$activos, bins = bins)
leaflet() %>%
addTiles(group="OSM") %>%
addProviderTiles(providers$Esri.NatGeoWorldMap , group="NatGeo") %>%
addProviderTiles(providers$CartoDB.DarkMatter, group="CartoDB-Black") %>%
addRasterImage(
fallecidos_rl,
opacity = 0.7,
group= "Fallecidos",
) %>%
addPolygons(
data = casos,
color = "black",
fillColor = ~pal(activos),
fillOpacity = 0.7,
weight = 1,
opacity = 1,
stroke = TRUE,
group= "Casos Activos"
) %>%
addLegend(pal = pal, values = casos$activos, opacity = 0.7, title = "Casos Activos") %>%
addLayersControl(
baseGroups = c("OSM", "NatGeo", "CartoDB-Black"),
overlayGroups = c("Fallecidos", "Casos Activos"),
options = layersControlOptions(collapsed = FALSE)
)